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We derive the non-abelian generalization of the Furry ap- 
proximation which describes the transverse dynamical evolu- 
tion of a hard projectile parton inside a spatially extended 
colour target field. This provides a unified starting point for 
the target rest frame description of the nuclear dependence 
of a large class of observables. For the case of the virtual 
7* — > qq photoabsorption cross section, we investigate then 
in detail under which conditions the nuclear dependence en- 
coded in the Furry wavefunctions can be parametrized by a q q 
QCD dipole cross section. The important condition is colour 
triviality, i.e., the property that for arbitrary AT- fold rescat- 
tering contributions the only non-vanishing colour trace is 
N c Cg. We give proofs for the colour triviality of the inelastic, 
diffractive and total photoabsorption cross section measured 
inclusively or with one jet resolved in the final state. Also, 
we list examples for which colour interference effects remain. 
Colour triviality allows us to write the 7* — > q q contribution 
to the DIS nuclear structure function F2 for small Bjorken 
XBj in terms of a path integral which describes the transverse 
size evolution of the qq pair in the nuclear colour field. This 
expression reduces in an opacity expansion to the N = 1 
result of Nikolaev and Zakharov, and in the eikonal approxi- 
mation to the Glauber-type rescattering formulas first derived 
by Mueller. In the harmonic oscillator approximation of the 
path integral, we quantify deviations from the eikonal limit. 
Their onset is characterized by the scales L/lf and E 1 " 1 L 
which relate the longitudinal extension L of the nuclear tar- 
get to the coherence length If and the total transverse energy 
E^ accumulated by the g-q-pair. 



I. INTRODUCTION 

The partonic interpretation of physical processes and 
their nuclear dependence is Lorcntz frame dependent. 
Drell-Yan, e.g., is a q-q —> 7*-fusion in the infinite mo- 
mentum frame [ |], but it becomes a 7*-bremsstrahlung 
radiation off the projectile quark if viewed in the tar- 
get rest frame [ ||, ||. Similarly, nuclear shadowing is 
due to the recombination of partons from different nu- 
cleons, if viewed in the infinite momentum frame [ Q. 
In the target frame, it arises from non-additive contribu- 
tions of the rescattering of the hadronic 7*-Fock states 
in a spatially extended medium [ ||, ^j, 0]. The de- 
scriptions in different Lorentz frames are equivalent, of 
course. But depending on the physical problem at hand, 
a well-chosen Lorentz frame may provide a particularly 
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simple partonic interpretation. To study the nuclear de- 
pendence of physical processes, the target rest frame pro- 
vides arguably the most intuitive picture. QCD factor- 
ization theorems can be expected to hold for the nuclear 
dependence of high-pt processes only [ |l0[ |ll|, [l^], and 
thus a large variety of other approaches to the nuclear 
dependence of hard processes exists in the literature [ 
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These exploit that in the target rest frame, the nuclear 
dependence of a physical process P with transition ampli- 
tude (\Jjjl f) can be attributed often to the multiple 
rescattering of the hard in- and outgoing partons inside 
the soft spatially extended target colour field. In the sim- 
plest cases, this rescattering effect can be described by 
multiplying the free in- and outgoing wavefunctions with 
straight eikonal Wilson lines which account for the lead- 
ing medium-induced colour rotation of projectile quarks [ 

In the present work, we derive and study a more re- 
fined approximation scheme for in- and outgoing wave- 
functions and \& / which includes the leading trans- 
verse dynamical evolution of the projectile partons in- 
side the target. To this end, we derive in section || ex- 
plicit expressions for and which approximate to 
leading order 1/E in the norm and next to leading or- 
der in the phase the solution of the Dirac equation in 
the presence of a spatially extended colour field. The 
1 /^-corrections included in these solutions are known to 
provide the leading contribution for observables which 
are essentially determined by the destructive interfer- 
ence between different production amplitudes as, e.g., the 
non-abelian Landau-Pomeranchuk-Migdal (LPM) effect [ 
[l| 0, [[! [tf}, 0. HI, H [H H or the nuclear depen- 
dence of Drell-Yan pair production [ ||, |TJ. Indeed, our 
derivation of \&j and 'F/ draws on exactly the same ap- 
proximation schemes which were used in recent studies [ 
p"3| , |27| of the nuclear dependence of these observables. 
In section ||, however, we discuss this rescattering effect 
without reference to a particular observable. In this way, 
we obtain a compact expression which turns out to be the 
non-abelian generalization of the abelian Furry approxi- 
mation [ |2^] and which may serve as building block 
in the calculation of very different nuclear dependencies. 

The main technical complication in working with non- 
abelian Furry wavefunctions is that they involve a path- 
integral over a path-ordered non-abelian Wilson line. In 



section III and IV, we show for the example of the 
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7* — > q q photodissociation process how explicit and ex- 
act calculations can be done despite this complication. 
The key step is a diagrammatic technique first used 



by Mueller and collaborators 
which allows to establish the colour triviality of certain 
cross sections involving TV-fold rescattering processes. 
Here, colour triviality means that the contribution to a 
medium-dependent observable to N-th order in the opac- 
ity involves only colour traces which reduce to the iV-th 
power of the Casimir. This renders the problem essen- 
tially abelian. In section [VC, we derive the correspond- 
ing diagrammatic identities for our configuration space 
formulation of partonic rescattering. We then use these 
identities in explicit proofs of the colour triviality of the 
inelastic and diffractive part of the inclusive and one-fold 
differential (i.e. one jet resolved) photoabsorption cross 
section. At least for the 7* — > q q photodissociation cross 
section, all colour trivial observables turn out to be given 
in terms of the transverse dynamical evolution of a QCD 
dipole which is described by a simple path integral /C. We 
discuss the eikonal limit in which the transverse dynam- 
ical evolution is neglected, and we quantify the leading 
corrections to this limiting case. 

We choose for all explicit calculations the photodissoci- 
ation process 7* — > qq mainly since it is the simplest ex- 
ample which allows to illustrate the main technical diffi- 
culties associated with the non-abelian Furry approxima- 
tion and colour interference effects. The 7* — > q q process 
is a quantitatively significant contribution to the nuclear 
structure function F2 , but phenomenological applications 
require the inclusion of processes with initial gluon radi- 
ation (e.g. 7* — > qqg). These are known to contribute 
in the aligned jet region, where the transverse qq sep- 
aration is large, to leading order in 1/Q 2 . Discussion 
of these radiative contributions involves the quark-gluon 
vertex and complicates the analysis of colour interfer- 
ence effects significantly. This lies beyond the scope of a 
first illustration of the non-abelian Furry approximation, 
given here. We shall address the corresponding addi- 
tional technical problems in a subsequent work where we 
study colour interference in the non-abelian LPM-effect. 
In the present paper, we discuss only two aspects of the 
very general question to what extent our approach can be 
applied to other or more exclusive observables or to other 



models of the medium: First, we give in section IV B and 



IV E examples of (more exclusive) photoabsorption pro- 
cesses which are not colour trivial, thus indicating the 
limitations of the approach advocated here. Second, we 
argue in section [v| that the specific model ansatz for in- 
medium rescattering is not instrumental for our results. 
In the Conclusions, we summarize our main results and 
we shortly comment on further perspectives. 



II. ASYMPTOTIC WAVEFUNCTIONS FOR 
RESCATTERING PARTICLES 

We want to calculate the nuclear dependence of some 
physical process P with transition amplitude P\^f) 
by describing how the medium affects the propagation 
of the in- and outgoing wavefunctions ^>i, For the 
case of the photodissociation process 7* — * qq depicted 
in Fig. 0a this means that we write the corresponding 
transition amplitude in the form 

P|%) =ie J d*s¥ u + (s,pi)7° 

xe^ 7 "e ik - x ^ v (x,p 2 ) . (2.1) 

Here, the physical process P is determined by the photon- 
quark vertex iee ■ 7 e lkx . In this section, we derive explicit 
expressions for the medium-dependence of these asymp- 
totic wavefunctions irrespective of a particular process. 



Section II A explains the solution of the corresponding 



abelian problem, section II B establishes the non-abelian 
extension. In the remaining sections of this paper, we 
shall discuss then for the example of the photodissocia- 



tion (2.1) how explicit calculations can be done. 
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FIG. 1. (a) Contribution to the photodissociation process 
7* — > qq with N- and M-fold final state rescattering for the 
outgoing quark and antiquark. (b) Sum over rescattering dia- 
grams for an outgoing antiquark which leaves the production 
vertex with momentum pi and colour 01. To leading order 
in opacity, this resums all rescatt ering effects and defines the 
non-abelian Furry wavefunction ( 2.12 ). 



A. Abelian Furry approximation 



We consider a relativistic positron with momentum 
(i?2,P2) (E2 3> m). Prior to its detection, this positron 
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undergoes multiple small-angle scattering in a spatially 
extended medium, described e.g. by a collection U(x) of 
single scattering potentials ipi(x) — <p(x — x*) localized 
at spatial positions x», 



(2.2) 



The asymptotic electron wavefunction ^ is a solution to 
the Dirac equation in this spatially extended field: 



d_ 

dt 



t/(x) — m 70 + i a ■ V 



*(x,p 2 ) = 0. (2.3) 



Which approximation for ty(x,p2) keeps the leading 
medium-dependence ? An expansion of \& in powers of 
the coupling constant does not since it amounts to an 
expansion in powers of the single scattering potential (p 
rather than containing the leading order effect of U. In 
QED, the leading [/-dependence is kept in the Furry ap- 
proximation [ ^2| which is a high energy expansion 
of the solution of the Dirac equation. For the outgoing 
positron wavefunction, it reads 



_ AE 2 t-ip 2 z 



£> 2 F(x,p 2 )?;(p2) 



(2.4) 



which is exact to order 0{U/E) and 0(l/E 2 ). Here, L>i 
denotes a differential operator 



Di = l-i 



2Ei 
a = 70 7 



2Ei 
n ■ x 



Pi 



Pi 



(2.5) 



and the unit vector n specifies the longitudinal direction. 
The differential operator acts on the transverse wavefunc- 
tions F. For very late times, i.e., for far forward longitu- 
dinal distances xl, this wavefunction satisfies plane wave 
boundary conditions 



F 00 (-K ± ,x L ,p 2 ) = exp <-ip2 ■ xj 



P2 
I ~ X L 



2 P 2 



F(y±,VL,P2) = J d*L±_ G(y±,y L ;xx,x L \p 2 
x F 00 (x ± ,x L ,p 2 ) ■ 



(2.6) 



(2.7) 



Its evolution to finite longitudinal distances is determined 
by the retarded Green's function G whose path integral 
representation reads (for z' > z) 



G{r,z;r',z'\p) = 

^2 



J Vr{0 exp | J d£ 



(0-<^(r(0,e) 



(2. 



Here, r = dr/d£ and G satisfies the boundary conditions 
r(z) = r, y(z') — r', with G(r, z; r', z' — z\p) — S(r' — r). 



The Green's function ( |2.8| ) describes the Brownian mo- 
tion of the projectile particle in the plane transverse to 
the beam. This is the leading medium effect on the prop- 
agation of the proj ectile. Starting from the abclian Furry 
wavefunction (2.4), the KST- formalism [ [H| ^2j then 
allows to describe e.g. the medium dependence of the 
LPM-bremsstrahlung spectrum [ p]| , p2[ . One of the main 
motivations for what follows is the question to what ex- 
tent the application of the same approach to non-abelian 
problems is justified. 



B. Non-abelian Furry approximation 

To find the non-abelian generalization of the Furry 
approximation, we consider a spatially extended static 
colour potential of the form 



i=l 



(2.9) 
(2.10) 



where T a , (a = 1, . . . , AT? — 1), denote the generators of 
the SU (N c ) colour representation of the projectile par- 
ton. The z-th scattering center is located at Xj and ex- 
changes a specific colour charge a = a,. 

For the spatial support of the potentials tpf we take 
the rapid fall-off of a Yukawa potential with some Debye 
screening mass M. We assume that the mean free path 
of the projectile parton in the medium is taken to be 
much larger than 1/M. This ansatz for ( 2.1 0| ) is known 
as Gyulassy-Wang model [ |l3|, Q and was originally in- 
troduced to mimic rcscattcring effects of hard partons in 
the colour-deconfined matter created in the early phase 
of a relativistic heavy ion collision. 

To leading order 0(1/ E), the Feynman diagram in 
Fig. |l|b is the only rescattering term. Rescattering 
contributions involving 3- or 4-gluon vertices are known 
to come with spatial suppression factors in the Gyulassy- 
Wang model [13, 17 . The JV- fold rescattering of a hard 



parton in the colour field (2.S) is thus determined by 



d 3 Pi ,3 *(A + m )7o 

(I X ? ; 



) 3 1 Pi ~ m2 + ^ e 
x [-iA (x i )]e- ix *-( p '+ I - p *)) «W(p). (2.11) 

Here, the quark leaves some production vertex with mo- 
mentum pi and colour a±. It undergoes N gluon ex- 
changes with the spatially extended colour potential A^, 
emerging in the "final state" with momentum p and 
colour a. I^ N \y) is a matrix in the colour representa- 
tion of the projectile. The path-ordering V implies that 
Ao(Xj+i) stands to the right of j4o( x «)- The diagram 
in Fig. ^ denotes the component la^l of this matrix. 
The momentum transfers to the quark line are written 
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as Fourier transforms of the static scattering potential 
with respect to the relative momenta p^+i — p^. Finally, 
instead of an e xplicit production vertex, we have intro- 
duced in ( 2.11 ) the incoming plane wave exp [— ipi ■ y]. 
This factor allows to glue the above Feynman diagram 
via y-integration onto another subprocesses P without 
specifying P at the present stage. 

In appendix [X[ we approximate the norm of i"W to 
leading order 0(1/ E) and the phase to next to leading 
order. Summing then over contributions for arbitrarily 
many N rescatterings in the medium, we find the non- 
abelian Furry wavefunction 

N=0 

= e tEy °- zpyL DF(y,p) w (r) (p) . (2.12) 
This expression compares directly to the abelian Furry 



approximation (2.4). In (2.12), D is the operator (2.5) 
with pi — p, and F is the outgoing transverse wavefunc- 
tion evolved from its asymptotic plane wave form to y 
with the Green's function G, 

F(y,p) = J dx±G(y±,y L ;yL±,x L \p) 

xf ffi (xi,i L ,p). (2.13) 

The Green's function G is the non-abelian generalization 
of the Green's function (2.8). It can be defined explicitly 
by its expansion in powers of the scattering potential Aq: 



G(r, z; r , z'\p) = G (r, z\ r , z\p) - i J d£ 

Z 

x J dpG (r,z;p^\p) A (p, £) G (p,^r',z'\p) 

XL XL 

+F J d£i J d& J dp 1 dp 2 G {r,z;p 1 ^ 1 \p) 

x i A (px , & ) G (p x , Ci; p 2 > 6 \p) 
xiA (p 2 ,&)G(p 2 ,fr,r',z'\p). (2.14) 

Here, Go is the free non-interacting Green's function 



G (r, z;t',z'\p) 



P 



a2 



2iri(z' — z) 



exp ■ 



ip (r — r 

%z> - z) 



(2.15) 



and the path ordering V in (|t|) ensures that the po- 
tential Aq(p 2 ,£,2) stands to the right of the pot entia l 
Aq(p 1 ,£ i i). We shall use the Green's function G in ( 2.14 ) 
for z' > z and z > z' by defining 

G(r', z'\ r, z\p) = G^(r, z; r', z'\p) , for z' > z . (2.16) 



For a hermitian scattering potential, this definition is 
compatible with the representation (2.14). Hermitian 
conjugation automatically inverts the path-ordering. A 



very compact notation for the expansion ( |2.14 ) can be 
given in terms of a path-ordered Wilson line VF([r]; z, z') 
which follows the non-abelian potential Ao from initial 
position (r(z),z) to final position (r(z'),z r ) along the 
path r(£) 

G(r,z;r',z» = 

= J 2>(£)exp j| |^r 2 (o| W([r];z,z'), (2.17) 

W([r]; Z) z') =PexpJ-i j d^M'tilt) j ■ ( 2 -18) 

Again, operators at larger longitudinal distances (z' > z) 
stand to the right. Expanding the Wilson line to fixed or- 
der O(Aq) coincides with the representation (2.14). The 
non-abelian Furry approximation thus differs from the 
abelian one essentially by path-ordering in the Green's 
function G which describes the dynamical evolution of 
the transverse wavefunction. 



To obtain the result ( 2.12 ), it is crucial that we approx- 
imate first the iV-fold rescattering diagram l( N *> , keeping 
the phase factor to order 0(1/ E) and then resum contri- 
butions for arbitrary N . If one parallels the derivation of 
the abelian case by directly starting from the solution of 
the non-abelian Dirac equation in the presence of a spa- 
tially extended colour field, one arrives at the recursive 
solution of Buchmiiller and Hebecker [ . The starting 
point for their recursion formula is a straight eikonal Wil- 
son line, and the recursion includes corrections to fixed 
order in energy. An explicit n-fold iterated expression 
is thus correct to fixed order 0(1/ E n ), but it does not 
contain the 0(l/i?)-corrections to the phase of the wave- 
function, which is characteristic for the Furry approxi- 
mation. Our diagrammatic approach allows to keep this 
phase. For large but finite incident energy, our solution is 
thus characteristically different from that given in Ref. [ 
[We have checked, e.g., that the abelian version of 
[ |29| does not lead to the correct medium dependence 
of the QED brcmsstrahlung spectrum while the abelian 
version of (2.12) does. This difference of both solutions 
is expected since the bremsstrahlung spectrum is sensi- 
tive to interference effects which stem from the 0(1/ E)- 
contributions to the phase of the wavefunc tions [ 

In the infinite energy limit, the solution ( 2.12 ) reduces 
to the straight eikonal Wilson line 



\im G(r,z;x ±1 x L \is) = W([r s }; z,x L ) , (2.19) 

(2.20) 



r a (f) = x ± hr 



x L - z x L - z 
which is the leading order input for the recursion formula 
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of Ref. [ 29 . This expression is often used to describe the 



leading colour rotation of the hard parton in a soft colour 
field [ HI M. 



The Furry wavefunction ( 2.12 ) derived here allows for 
the description of the nuclear dependence of a large class 
of physical observables. To discuss its limitations, we 
note that QCD factorization theorems for observables 
including rescattering effects have been established (and 
can be expected to hold) only for processes for which 
kinematical constraints ensure very high transverse mo- 
mentum transfers (> A QC d), [ |§. The de- 
scription of soft nuclear rescattering of hard partons al- 
ways depends on additional assumptions and different 
approaches may be taken [ |l|, [l3| |jq| . In particular, we 
assume (as all other treatments do) that the leading nu- 
clear dependence can be obtained using the idealization 
of "asymptotic" parton wavefunctions, i.e., without con- 
sidering parton fragmentation. Since parton fragmenta- 
tion and soft parton rescattering involve transverse mo- 
menta of the same order, this is a requirement on the 
spatial separation of the two phenomena and amounts to 
an energy-dependent limit on the longitudinal extension 
L of the nuclear target up to which this description can 
be expected to apply. Moreover, the ad hoc separation 
of the transition amplitude into a production P and a fi- 
nal state wavefunction including rescattering effects may 
be oversimplified for some processes. We think, e.g., of 
the production and rescattering of a heavy quarkonium 
state with specific quantum numbers. Taking P to be 
the production of the heavy q-q-paiv, it is not clear a pri- 
ori to what extent final state rescattering modifies the 
quantum numbers of this state. First studies indicate [ 
that the discussion of this problem requires a classi- 
fication of the hardness of medium-induced momentum 
transfers which lies outside the scope of the present cal- 
culation. To sum up: the discussion of the nuclear depen- 
dence of hard observables in terms of (Q^lPl^//) seems 
justified if no other soft scale (introduced e.g. by par- 
ton fragmentation or by the binding energy of the final 
state) interfers significantly with the final state rescat- 
tering effect described by the Furry wavefunctions. The 
non-abelian LPM-effect [ [u| 0, pL [If] , the nu- 
clear dependence of Drell-Yan yields [M El] and nuclear 
shadowing [ ^, [?| are prominent examples for which a 
description in terms of (^^P^f) seems suitable. 



III. PHOTODISSOCIATION 

In this section, we give a closed expression for the 
7* — > qq photodissociation cross section cr 7 _>9 ' in terms 
of the non-abelian Green's function G and the squared in- 
going wavefunction $ of a freely evolving q-q-p&vc. This 
photodissociation cross section contributes to the total 
virtual photoabsorption cross section o-J otal which is re- 
lated to the deep inelastic structure function F2 , [ |3l| 



F 2 (x,Q 2 



Q 2 



47r 2 a c 



tal ( X > 



(3.1) 



In the target rest frame, one can show that photodissoci- 
ation as depicted in Fig. ||b dominates for small Bjorken 
x over the j*-q fusion process shown in Fig. ||a [ 
The leading corrections to c^ ot ^ 99 come from initial state 
gluon radiation [ fj], [33| which we neglect in what follows. 

We consider a virtual photon of four momentum 
k = (is, 0j_, yfv 2 + Q 2 ) which dissociates into a quark- 



antiquark pair in the time interval 
The dissociation cross section is 

-v*-^™ 1 



-T/2 < t < T/2. 



'total 



dX\S 



2\k L \T' "- v ■ (3 - 2) 

where the phase space element is written in terms of the 
on-shell momenta of the outgoing quarks, 

_ (2^)3 2)^1 (27r)3 2|^|- ^ 

The rescattering of these quarks in the spatially extended 
colour field is described by the Furry wavefunctions in the 
dissociation amplitude 

Sfi = S%e», (3.4) 



c</< 



xe- e ^e ik - x ^ v (x,p 2 ) 



(3.5) 



To discuss the case of different photon polarizations, we 
write this disscociation amplitude as a contraction of 
with the polarization vector e M . To shorten our notation, 
we do not keep track of the fractional charge e q of the 
quark and of the number of flavours. To do this, all our 
final results have to be supplemented by a sum over the 
available flavour channels weighted by e 2 . 



(a) 




time 



(b) Photon 




FIG. 2. Two possible time orderings for the interaction of 
a virtual photon with a nuclear target. At small Bjorken x, 
the 7* — > qq photodissociation process (b) dominates over the 
7*-g- fusion, see Ref. [ m. 



5 



Inserting the Furry wavefunctions for quark and an- 



tiquark, see (2.12), we can separate in ( |3 . 5[ ) an energy 
conserving 8- function, 5£ = ie (2n) 8{E 2 — E\ — u) M^. 
We introduce the fraction a of the initial energy carried 
by one quark, E2 = a v, E± = — (1— a) v. Here, E\ is neg- 
ative, since the corresponding 4- momentum p\ flows into 
the quark-photon vertex, see Fig. [j]a. After coordinate 
transformation (pf ,p 2 ) — > (a,k L ), the photodissociation 
cross section reads 



' total 



da 



Av 2 a (1 - a) 
dpi dpi 



(2tt) 2 (2tt) 2 



(\M ft \ 2 ) 



(3.6) 



Here, we have used Fermi's golden rule and the k L - 
integration to eliminate the two energy-conserving 8- 
functions. The brackets (...) denote an in-medium aver- 
age over the colour field which will be specified below. 
The probability |M/i| 2 can be written in terms of the full 
interacting Green's functions (2.17) as 



\M fi \ 2 = / d 3 yd 3 y e i 9(t/X,-w) e -e(|!/z,| + |S i |) / dx± 

xrfx'_L dx ± dx' ± e *Pi<*±-*±) e-iP'-^-*^ 
xG(x;y\p 2 )e^G(y;x'\ Pl ) 

xG(x';y\p 1 )e„T* U G(y;x\p 2 ) . (3.7) 

In general, all spatial coordinates can be different for 
Mfi and M'^. We characterize those for Mj i by a bar. 
The spinor st ruct ure of the Furry wavefunctions ( p. 12 ) is 
contained in ([Tt]) in the vertex function 



f " - u^\- Pl ) b\ 7 ° y b 2 «W(p 2 ) 



(3.8) 



In appendix [b], we discuss how these vertex functions 
combine with free Green's functions Go to the incoming 
q-q Fock state. In terms of the square $(Az; Az; a) of 
this Fock state, the differential photodissociation cross 



section (3.6) takes the form 



total 



da dpi dp^ 

= (S* / dhl dh2 ^ dS 2 *(Ab; Ab; a) 
x J dx ± dx' ± dx ± dx' ± e 4p i L - (x ' J 
x(G(b 2 ;x| J52 )G(x;b 2 |p 2 ) 
xG(b 1 ;x , |p 1 )G(x';b 1 |p 1 ) 



-ip^--(xj_-X_L 



(3.9) 



where Ab = bi b 2 and Ab = bi b 2 . The explicit 
form of the squared incoming wavefunction $ is derived 
in appendix M. 



In equation ( |3.9| ), we assume that the nuclear target 
has nonvanishing density only for some longitudinal po- 
sitions zl > 0. The integration variables b i; bj denote 
the transverse boundary values of the Green's functions 
at the front end Zl = of the nuclear target. Fig. || 
gives a graphical representation for ( |3.9| ) which we shall 
heavily use in what follows. According to Fig. ||, the 
virtual photon wavefunction starts interacting with the 
nuclear medium at zl = in both the amplitude and 
complex conjugate amplitude after dissociating at longi- 
tudinal positions j/l, yi, with uliVl < 0. We em phasize, 
however, that despite appearance, equation (3.9) as well 
as its graphical representation in Fig. || also include the 
case that the photon vertex dissociates inside the target 
at ul > 0. Technical details of how this is handled are 
given in appendix H. 



>'L 



CUT 




FIG. 3. Re pres en tatio n of the 7* — > qq photodissociation 
cross section (3.6), (3.7). In the amplitude Mfi, the photon 



dissociates at position (y,J/i) into a quark and an antiquark. 
These propagate with free non-interacting Green's functions 
(dotted lines) up to the front end of the target at longitudinal 
position zx. From zl onwards, the quark and antiquark are 
propagated inside the colour field with interacting Green's 
functions (full lines) up to their asympotic on-shell states 
("cut"). For the complex conjugate amplitude, the analo- 
gous evolution of the hadronic qq Fock state is depicted on 
the other side of the cut. 



IV. OPACITY EXPANSION AND COLOUR 
TRIVIALITY 

In this section, we analyze the photoabsorption cross 
section (3.9) by expanding the interacting Green's func- 
tions G in powers of the scattering potential ^4o( x )- In 
a first step, we focus on the total photoabsorption cross 
section 



°-lt7i qQ = "cm / da I dbi db 2 dbi dh 2 $(Ab; Ab; a) 



dx±_ G(b 2 ;x|p 2 ) G(x;b 2 |p 2 ) 
• / dx' ± G(b 1 ;x'b 1 )G(x';b 1 b 1 )). (4.1) 
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To make sense of an expansion in powers of Aq, we have 
to specify the in medium average (. . .) in terms of Aq. 



To this end, we specify for the colour potential (2.9) the 
contribution of a single scattering potential, centered at 



tf(x J . ) = « oo < 



d 3 K 

(2tt) 3 



(4.2) 



Here, we have approximated the argument of the single 
scattering potential ao(/«) « o,q(k±). This implies that 
the momentum transfer occurs at a fixed longitudinal 
position £ = £j. It is the standard approximation in the 
high energy limit where the dominant momentum trans- 
fer from the medium is transverse, and it motivates the 
use of time-ordered perturbation theory to rescattering 
problems [ fD| pp. 

Starting from ( |4.2| ), we define the medium average (. . .) 
as an average over the transverse and longitudinal posi- 
tions (fj,ij) of the scattering potentials and the colour 
factors a*. 



x/(fi, ...,i N ;zi,.. .,z N ;ai, . . .,a N ) . (4.3) 

Here, A± is a total transverse area which we divide out 
to regain the cross section per unit transverse area. N is 
the number of different single scattering potentials up to 
which the function / is expanded. 

To simplify notation, we replace in what follows the 
discrete sum over Zj in ( |2.9| ) by an integral over the den- 
sity n of scattering centers 



A (x ± ,£)= dz l n(z i )<p?(x±,OT a . 



(4.4) 



Since the effective momentum transfer from the single 
potential ipf(x.±,£) occurs at £ = ij, we shall often work 
with £ as integration variable. The expansion of the 
photoabsorption cross section to 7V-th order in Aq is an 
expansion in the 7V-th order of the opacity parameter 
a 2 . J d^ n(£). In what follows, we study the perturbative 
expansion in this parameter. 



A. N=l result of Nikolaev and Zakharov 

As a first illustration of the above formalism, we ex- 
pand th e in tegrand of the total photoabsorption cross 
section (hi) to first order in the opacity parameter 
a 2 J d^ n(£), thereby reproducing the well-known nuclear 
shadowing result of Nikolaev and Zakharov. For N = 1, 



The zeroth and first order contributions O(Aq) and 
O(Aq) to crl ^d 9q vanishes due to energy-momentum con- 
servation: without momentum transfer to the medium, 
the q-q-pSLiT cannot appear on-shcll. To second order 
O(Aq), there are four different terms which we depict 
diagrammatically in Fig. |Ja: for each term, the full lines 
correspond to the full lines shown in Fig. |^ and describe 
the full Green's functions in the total photoabsorption 
cross section ( |4.l| ). To simplify the representation, we 
have dropped in comparison to Fig. ^ the photon lines 
and incoming q-q wavefunctions. For the total photoab- 
sorption cross section, the transverse positions at the cut 
are equal for the amplitude and the complex conjugate 
amplitude. 




FIG. 4. (a) Diagrammatic representation of the N — 1 
contribution to the total photoabsorption cross section ( |4.l[ ). 
For each term, the upper (lower) line denotes the Green's 
function with energy p\ (p?) and the cusps denote the position 
of the cut. (b) Definition of th e diagrammatic shorthand used 
in (a) and calculated in (4.5). Crosses stand for one power 
of the scattering potential Aq. The medium average ensures 
that the momentum flow k through the cut is conserved and 
that the potential Aq is linked to the q-g-system in Ma and 
Mff at the same longitudinal position. 

Let us consider the first term on the r.h.s. of Fig. |^. 
Here, the Green's functions of argument p 2 are free while 
both Green's functions with argument p\ are expanded 
to first order. Averaging over the position (fi , Z\) of the 
center of the scattering potential, we find 



we have to expand all Green's functions in (4.1 ) and then 
to collect all contributions to order 0(Aq). 



dxj_ G (b 2 ;x|p 2 ) G (x;b 2 |p 2 ) 

x J dx , ± GW(bi;x / |pi)Gf 1 )(x / ;bi|piJ 

J <2£in(£i) J -^^\ao(K ± )\ 2 T ai T ai 
x,5( 2 )(b 2 -b 2 )(5( 2 )(b 1 -b 1 ). (4.5) 



This term does not depend on the separation between 
the q-q pair since the potential touches only one of the 
quark lines. In contrast, the second term on the r.h.s. of 
Fig. ||a depends on this separation and takes the explicit 
form 
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dem(6) 



(K ± )| 2 T Ql T ai e-^" (bl - b2) 



xe 



^fe>^--) ^ 2 > (b 2 - b2 - - 

(4.6) 



x<5< 2 > (bi-bi-^^-^)) 

This shows explicitly that the transverse Brownian mo- 
tion induced by rescattering is taken into account in our 
formulation: the q-q dipole does not propagate at fixed 
transverse separation bi — bi but changes its size in re- 
sponse to an interaction at the longitudinal position £1 by 
terms of the form ^-(fi — Zl)- For N = 1, the first and 
only interaction ta&s place by definition at the position 
zl = £1 and these terms vanish. For N > f, however, 
they do not vanish and lead to a nontrivial dynamical 
evolution of the transverse dipole size, see below. 

Adding up the four real contributions in Fig. |]a, and 
doing the colour trace, we find 



total 



N c a e 



t^n(0 / da 



x J db$(b;b;a) a(b) , (4.7) 



a(b) = 2C F 



dn± 



\a (K ± )\ 2 (l-e-^' b ) . (4.8) 



This is the result of Nikolaev and Zakharov [ |) . The total 



cross section (4.1) is the product of the Born probability 
<£>(b;b;a) for the incoming hadronic Fock state times a 
dipole cross section <r(b), integrated over the transverse 
size b of the dipole and the energy distribution a between 
the quark and antiquark. 



B. N = 2 tagged: non-trivial colour interference 
effects 

For more than N = 1 scattering center, there are in 
general non- vanishing interference terms between ampli- 



tudes containing different powers of Aq, see section IV C 



below. Such contributions appear in the generic situa- 
tion in which one knows the distribution of scattering 
centers but no information is available on whether these 
centers have participated with finite momentum transfer 
in the scattering process. At least in a gedanken exper- 
iment, however, this additional information can be ob- 
tained from the recoil of each scattering center, i.e., by 
tagging the scattering centers. For this very exclusive ob- 
servable, all contributions to the cross section contain the 
same power of A in both amplitudes, and all scattering 
centers contribute to the amplitude with finite momen- 
tum transfers. We say, such scattering centers contribute 
with real terms, in contrast to the contact term contribu- 



were considered recently in the study of the non-abelian 
LPM-effect [ Q , and there is some interest in comparing 
results of tagged and untagged scenarios to understand 
the origin of colour triviality. To this aim, we consider 
here as a second application of the opacity expansion the 
case of two tagged scattering centers: (i) two scattering 
centers are involved in the process (4.1) and (ii) both 
scattering centers provide real terms, i.e., transfer finite 
transverse momentum on the amplitude level. The 16 
contributions to this ^7agged'(-^ = 2) are listed in Fig. ||. 
Summing them up, we find 

dem(fi) J d£ 2 n(£ 2 ) J da J ' db 

x (Tr[aabb}A (1) +Tr[abab] A (2) ) , (4.9) 

where we have introduced the notational shorthands 

= 4$(b,b;cv) (l - e - tKl ± b ) , (4.10) 

A^ = 4 * f b - K2 ' ±( f ~ M , b + " 2 ^ 2 ~ & > ; g) 
\ 2va(\-a) ' 2i/a(l-a) / 



3 — I Kl ,_L'/«2, 



-»2,j 



(4.11) 



tions which we discuss in section [V C . Tagged scenarios 



Here, the term A^ contains information about the non- 
trivial dynamical evolution of the g-g-separation which 
is determined by the "relative transverse mass" /i = 

— — ) = va(\ — a). The momentum transfer Ko i 

Pi P2 J y ' Z '- L 

determines this change of the g-g-size via a classical equa- 
tion of motion: oc K2 ^ ± (£2 — £i)- I n the high energy limit 
v — ► 00, this transverse motion can be neglected and the 
g-g-pair propagates at fixed separation. We find 

lim ol^{N = 2) 

= a om J d£m(£i) J cZ£ 2 n(£ 2 ) J da J db 

f dK% I dK 2 ,± , , \ 1 2 1 / m2 

X J "(2^) T "(2^ |ao(Kl ^ )l |ao(/t2 ^ )l 

x$(b,b;a)4 (l - e - iKl > x,b ) 

x (Tr[aabb}-Tr[abab}e- lK - 2±h ) . (4.12) 

Even in this limiting case, the colour algebra does not fac- 
torize from the momentum dependence. This is a con- 
sequence of the non-trivial colour interference between 
different photodissociation amplitudes. The photoab- 
sorption cross section cannot be written as a function 
of the dipole cross section a(b). We note that this com- 
plication stems entirely from the non-abelianess of the 



s 



problem. If we replace both colour traces in (4.12) by 
the same constant, the two /^^-integrations combine to 
the square of the dipole cross section ( |4.8| ). To sum up: 
in the tagged case, non-trivial colour interference effects 
prevent us from representing the photoabsorption cross 
section as a function of the dipole cross section. 







FIG. 5. The 16 contributions to the N = 2 photoabsorp- 
tion cross section (4.£). In this "tagged" case, exactly two 
real momentum transfers occur on the amplitude level. The 
tagging gives rise to non-trivial colour interference effects. 



C. Contact terms 



In the calculation of the N = 1 cross section (4.7), we 
did not include all terms of the Taylor expansion of (4.1) 
to order O(Aq). A second order contribution which was 
not included is e.g. 



dx± G (b 2 ; x\p 2 ) G (x; b 2 \p 2 ) 

x j dx / ± G (b 1 ;x / |p 1 )G( 2 )(x';b 1 b 1 )J 

-~ jd^n(^) J ^\a ( K± )\ 2 T ai T at 

xS^(b 2 -b 2 )S^(h 1 -b 1 ). (4.13) 



Here, the potential Ao is linked two times to the com- 
plex conjugate amplitude and zero times to the ampli- 
tude. The medium average (. . .) results in two impor- 
tant constraints: (i) both powers of Aq couple to the q-q- 
pair at the same longitudinal position £i, and (ii) no net 
momentum is transferred to the q-q- system. Scattering 
contributions with these properties were called virtual by 
Mueller and collaborators [ . We refer to them as 

contact terms, all other interactions are refered to as real. 
Diagrammatically, the six contact terms for the N = 1 
case are given in Fig. ^a, and equation (4.13) is repre- 
sented by the first diagram on the r.h.s. 







FIG. 6. (a) The six contact terms for the N = 1 scattering 
cross section, (b) Definition of the diagrammatic shorthand 
used in (a). The medium average ensures that (i) the momen- 
tum flow k is conserved and no net momentum is transferred 
to the g-q-system, and that (ii) both powers of the scattering 
potential act at the same longitudinal position z* . 

For the N — 1 photoabsorption cross section ( J4. 7| ) , con- 
tact terms do not contribute since at least one interaction 
with the g-g-system is needed in both the amplitude and 
complex conjugate amplitude in order to get the final 
state on shell. For N > 2, the same argument about 
energy-momentum conservation does not imply the ab- 
sence of contact terms. We have to distinguish the fol- 
lowing cases: 
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'inol. 

one real interaction. 



1. CjncT"^ inelastic: the cross section involves at least 



2- <4 ff 



qq diffractive (or elastic): the medium inter- 
acts with the g-g-pair only via contact terms. 



3. cr^^ 99 total: the cross section involves at least one 
(real or contact) interaction in both Mfi and M^. 

The notion "inelastic" is justified since any real inter- 
action changes the colour of the target and thus affects 
the hadronic activity between projectile and target rapid- 
ity (in this sense, the target "breaks up"). The notion 
"diffractive" or "elastic" involves an additional assump- 
tion: by construction, contact terms transfer exactly zero 
net transverse momentum to the g-q-pair. Also, we re- 
gard in all calculations and for arbitrary scatterings the 
longitudinal momentum transfer as sufficiently small to 
justify the approximation ao(«) « ao(n±). Our assump- 
tion in the above definitions of the diffractive and total 
photoabsorption cross sections is that despite neglecting 
the longitudinal momentum transfer in all explicit cal- 
culations, it is sufficient to put the (extremely weakly 
off-shell) partonic contributions of 7* on-shell. This is 
usually assumed in the calculation of the total [ [j| and 
diffractive [ [| cross sections. 

For N > 1, contact terms obviously play an important 
role in the calculation of the inelastic, diffractive and 
total photoabsorption cross section. In the remainder of 
this subsection, we summarize some of their important 



properties and shorthands, which will be heavily used in 
the following subsections: 



1. Identities involving contact terms 



factor Cf, i.e., we can substitute in (4.14) 



First we observe that in ( 4.13 ), the Green's func- 
tions of momentum p 2 result in the transverse 5-function 
S^(h 2 — b 2 ). Dro pping the m on both sides of ( 4.131 ), 
and observing that ( 4.13 ) and (4.5) differ by a factor 
we have checked the identity in Fig. [7|a. By a similar 
calculation, one also checks the identity Fig. 0b. 



(a) 

(b) 
(c) 
(d) 
























FIG. 7. At the cut, there are simple identities between real 
interaction terms and contact terms (a), and amongst real 
interaction terms (b). (c) A combination of the three possi- 
ble contact terms at the same longitudinal position £ leads to 
a dipole cross section a(ri(£) — r2(£)) whose trans verse size 
is evolved further with free Greens functions, see (4.14) (d) 



Iterating the combination of contact terms (c) leads to a dy- 
namical evolution of the dipole size between initial transverse 
size ri — r 2 and final size r j — r' 2 . This evolution is described 
by the path- integral (4.15). 



2. A path integral from iterating contact terms 

We now consider the case in Fig. 0c, where the q-q- 
state is evolved with exactly one contact term between 
longitudinal positions z and z' from a transverse separa- 
tion ri — to a transverse separation r[ — r' 2 . The sum 
over the three possible contact terms at longitudinal po- 
sition £, represented on the r.h.s. of Fig. 0c, takes the 
explicit form 



1. 



Then, equation fl4.14|) is seen to correspond to the first 
order n(£) density expansion of the double path integral 
(for z' > z) 

M(r l ,r 2 ,z;r' 1 ,r 2 ,z'\a\u) 

X>ri J Vr 2 

ri(z)=ri r 2 (z)=r 2 



x exp < i 



di-{avl + {l-a)v\) 



+il„(0^(r 1 (0-r a (0)}. (4.15) 



Moreover, the expansion of (4.15) to m-th order in the 
density n(£) corresponds exactly to the m-th term on the 
r.h.s. of Fig. 0d. This shows that the iteration Fig. 0d of 
the one contact term interaction Fig. 0d is described by 
the path- integral ( 4.15| ). Since this path integral plays 
an important role in what follows, we simplify it further. 
Using the coordinate transformation 



r Q (0 = (l-a)n(0+ar 2 (0, 
r 6 (0 = ri(0-r 2 (0, 



(4.16) 
(4.17) 



we can write 



M(ri,r 2 ,z;r' 1 ,r' 2 ,z'\a\u) 
= JCo(r a {z'),z';r a (z),z\v) 

xK,(r b (z'), z'\ r b (z), z\va[l - a)) 



(4.18) 



The path-integral on the r.h.s. of this expression is given 
by (for z' > z) 



K,(r(z'),z';r(z),z\fj) 
= / Vr exp < i d£ 



tlr 2 + i-n(0a(r) 



(4.19) 



Its opacity expansion reads 
K(r', z'; r, z) = /C (r', z'; r, z) 



T a ---T a 
2C f 



[d£n(£) y^OdratOo-MO-raCO) ~/ ^/ rf /°^o(r', z'; r, ^) ^) /C (p, ^; r, 



x G (r 2 , z; r 2 (£), CIP2) G (r 2 (O, 6 r 2 , z'\p 2 ) 

x G (ri, z'- n(0, i\vx) G (ri(0, 6 r 2) z\ Pl ) . (4.14) 

Since Fig. 0c is only part of a more complicated Feynman 
diagram, we can say nothing about the positioning of the 
colour factors T a . The corresponding expression T a ■ ■ ■ T a 
in (4.14) is purely formal. However, for the following, we 
anticipate an argument explained in the next subsection: 
for our purposes, the colour factors reduce to a Casimir 



+ / ' d£i J d& J d Pl dp 2 /C (r',z';p 2 ,£ 2 )£(p 2 ,£ 2 ) 

x/C(p 2 , 6; Pi, £1) S(pi, Hi) /C (pi, £1; r, z) , (4.20) 

where £(p, £J = ^n(£)<r(p). Here, we have suppressed 
the explicit /^-dependence in /C. The corresponding free 
Green's function /Co reads 
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K. (r',z';r,z) = 



2ni(z' - z) 



exp ■ 



ifi (r' 



2 0' - z) 



(4.21) 



In analogy to our definition of the Green's func- 
tion G in ( 2.16Q , we define K.{r(z), z; r(z'), z'\fij = 
K,Ur(z'),z';r(z),z\fx) for z' > z. 

A path integral K, of the form derived here was first 
used by Zakharov in the abelian problem of calculat- 
ing the passage of ultrarelativistic positronium through 
matter [ |3^]. More recently, the same path integral 
was shown to appear in the final expression of the 
medium-dependence of the LPM-bremsstrahlung spec- 
trum [§]| ii . 



D. N = arbitrary, untagged: contact terms remove 
colour interference effects 

What happens if contact terms are included in the 
N = 2 calculation of the elastic, inelastic and total pho- 
toabsorption cross section ? In Fig. ||, we have classified 
into four subsets all terms of order 0(n(£) 2 ) which have 
to be considered for the N — 2 photoabsorption cross 
sections. There are in Fig. ||a the 4 x 4 = 16 terms for 
which both interactions are real, in Fig. ||d the 4 x 6 = 24 
terms for which only the second interaction is contact, in 
Fig. ||c the 3x4x2 terms for which only the first inter- 
action is contact, and in Fig. ||d the 2 x 3 x 6 = 36 which 
involve two contact interactions. All together, these are 
100 terms. For the total and the inelastic cross section, 
the terms in Fig. ||a and Fig. ||b both contribute, but they 
cancel each other exactly. More precisely, the first dia- 
gram in Fig. ^a cancels the first diagram in Fig. ||b, the 
second diagram in Fig. ||a cancels the second diagram 
in Fig. ||b, etc. This is a consequence of the identities 
Fig. 0a and which we apply to the last interaction 
before the cut. 

The same argument cannot be made for the cancel- 
lation of the diagrams in Fig. He and Fig. ||d: energy 
momentum conservation restricts the occurence of dou- 
ble contact terms, since at least one interaction is needed 
on the amplitude level to get the final state on-shell. De- 
pending on whether we calculate the diffractive, inelastic 
or total photoabsorption cross section, different combina- 
tions of these two sets of diagrams contribute. The real 
interaction Fig. |^a at the cut is determined by minus 
the combination Fig. ||a of contact terms. This allows 
us to express all three contributions to the photoabsorp- 
tion cross section in terms of contact terms, as shown in 
Fig. H Most importantly, all contributions to these pho- 
toabsorption cross sections have the colour trace Tr [aabb] : 
for the untagged case for which contact terms are in- 
cluded, colour interference terms vanish in the diffractive, 
inelastic and total N — 2 photoabsorption cross section. 




b =0=0= 



(b) 



" =0= + =0= 



(d) 



(e) inelastic = (c) 






diffractive = (d) = 



total = (c) + (d) 



FIG. 8. The contributions to the N — 2 photoabsorption 
cross section: (a) the 16 real contributions which are present 
in the "tagged" case, (b) Contributions for which the first 
interaction is real and the second a contact term, (c) Contri- 
butions for which the first interaction is a contact term and 
the second interaction transfers momentum, (d) Terms which 
contain two contact interactions. Here, energy-momentum 
conservation ensures that the two contact terms do not stand 
on the same side of the cut. (e) Inelastic, diffractive and 
total photoabsorption cross section are determined by differ- 
ent combinations of the contributions (c) and (d). Rewriting 
the real interaction in (c) as the sum of two negative con- 
tact term contributions, a particularly simple representation 
is obtained. 

The above argument can be generalized easily to ar- 
bitrary N. As shown in Fig. |9|a, whenever at least one 
of the (N — 1) first interactions is real, the sum of all 
real and virtual terms for the N-th interaction cancels. 
Inelastic, diffractive and total photoabsorption cross sec- 
tion can thus again be described in terms of (N — 1) 
contact terms and a characteristic iV-th contribution at 
the cut. As an immediate consequence, the colour trace 
of all these contributions reads 



Trp^aj T a2 . . . T aN T aN T a ^ 



.T ai ]=N c C£. (4.22) 



In contrast to the tagged photoabsorption cross section, 
the colour structure factorizes from the momentum de- 
pendence: Inclusion of the contact terms results in colour 
triviality for the iV-fold diffractive, inelastic and total 
photoabsorption cross section. This justifies the assump- 
tion made in deriving the path integral (4.15) in the last 
subsection. 
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I 











+ 




(a) 



(b) inelastic = 



(c) diffractive = 



(d) total = (b) + (c) = 



FIG. 9. (a) Cancellation of real and contact terms for any 
combination of interactions [I] which contains at least one real 
or contact interaction, (b)-(d) The identity (a) allows for the 
generalization to arbitrary N of the N — 2 photoabsorption 
cross section of Fig. ^. 

According to Fig. ||(b)-(d), the inelastic, diffractive 
and total photoabsorption cross section can be written 
in terms of the path- integral JC of ( 4.19 ) which describes 
the iteration of contact terms. For the inelastic photoab- 
sorption cross section, this leads to the expression 





il<\ I dh dh 

L 



d£n(0 / dr(0 



x$ (b; b; a) 



x£(b,0;r(O,£Ml-«)Mr(O) 
x£(r(f),£;b,0|z/a(l -a)) . (4.23) 

Here, the Green's functions JC describe the propagation of 
the g-g-dipole from the front end of the target at z — 
up to the position £ of the last interaction. The lower 
boundary £ n i of the integration over £ is defined by the 
position of the next to last ("nl") interaction. This is an 
awkward property, since we have to expand the ZC's in the 
number N of interactions to determine £ n i ord er for order 
in N. This makes it important to find for ( 4.23 ) with 



the help of the opacity expansion ( 4.2C ) the equivalent 
representation 



^ineT 9q =N c a om Ida I dbdb$(b;b;a) 



<5< 2 > (b - b) 

dr o ^(b,0;r o ,i||i)^(r ei i;b 1 0|/<) 



(4.24) 



For the diffractive contribution, one finds 



Miff. 



N c a cm da / dh dh $ (b; b; a) / dr. 



S (2) (b-r e ) -£(b,0;r e ,L|^ 
5 (2) (r e -b)-/C(r e ,L;b,0|M) 



(4.25) 



The total photoabsorption cross section is the sum of 
both 



7 i _ N a 
a total — Jv c a c 



da / dbdb$ (b; b; a) 



2 5 {2) (b-b) -£(b,0;b,L|/*) 
-/C(b,i;b,0|/x) 



(4.26) 



Before turning in the next subsection to the discussion 
of these expressions, we remark shortly on a phase con- 
vention implicitly used in the above results: 

The first (N — 1) interactions in the above photoab- 
sorption cross sections are contact terms. Contact terms 
stand either to the right or to the left of the cut. The first 
interaction thus occurs for each non-vanishing contribu- 
tion to (4.24) at different longitudinal positions z a , Zb in 
the amplitude and complex conjugate amplitude. As a 
consequence, the free incoming wavefunction has to be 
evolved to these different positions in Ma and Mg. We 
explain in Appendix |b] following (B8) that the expression 
$ for the squared incoming wavefunction is only correct 
as long as z a = Zb- If one insists on using <i> for z a ^ Zb, 
one is forced to adopt the phase convention: 

J e?b $ (b; b; a) /Co(b, z a ; r, Zb\va(l — a)) 

= $(b;r;a) e - * 9 . (4.27) 



Strictly speaking, using $ one has done the i/l- and yi,- 
integrations in the photodissociation probability (j§7i| ) be- 
fore specifying the true endpoints of these integrations 
which are given by the positions z a , Zb of the first in- 
teractions in Mg and Mg. The phase convention ( 4.27 ) 
corrects for the part of the i/l integral which one misses 
in assuming z a = Here, we shortly illustrate the con- 



sequences of ( 4.27 ): 



With the help of (4.27), one checks immediately that 
the N = 1 contribution to (4.24) agrees with the result 



(4.7) of Nikolaev and Zakharov. There is no additional 
phase factor in this case since z a — zi,. However, there 
is an additional phase factor for all contributions N > 1. 
In the notationally simplest case, N = 2, we find 



ji -*gg 

inel. 



(TV = 2) 

-N c a cm Re I da I d? n(?) I d£n(£) 



drdr$(r;r;a) e 4<z (^') 
er(r) £q (r,f';r, £!/■«) 0"(r) . 



(4.28) 



The phase factor exp{iq — ^')} is sensitive to the differ- 
ence between the longitudinal position of the first scat- 
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tering center in Mg and Mg. The scale is set by the 
inverse coherence length 



91 

2v 



2va{\ 



1 



(4.29) 



It can be neglected in the high energy limit v — > oo where 
the longitudinal extension of the target becomes small 
compared to l/q. 



E. Colour triviality of the one-fold differential 
photoabsorption cross section 

Colour triviality of the inelastic, diffractive and total 



photoabsorption cross section (4.24)-(4.26) is the result 
of a complete diagrammatic cancellation which arises due 
to the identities in Fig. |?] (a) and (b) . These identities are 
stronger than implied by the optical theorem since they 
are based only on the transverse momentum integration 
of one of the two quarks in the final state. This makes it 
possible to use them in the simplification of more differen- 
tial observables. Here, we study the one-fold differential 
photoabsorption cross section 



da' 1 a c 



dhi db 2 dbi db 2 $(Ab; Ab; a) 



da dpi (2tt) 2 

x / dx JL dx'xdx'_ L e ip tix'±-*'jL') 



x^G(b 2 ;x|p 2 )G(x; b 2 \p 2 ] 
xG(b 1 ;x , |pi)(5(3E';bilpi; 



(4.30) 



We shall derive for the inelastic part of (4.30) an expres- 
sion in terms of 



cr(b) = 2 C F 



dn± 



\a (K±)\ 2 e 



-ki b 



(4.31) 



which is closely related to the dipole cross section, 

a(T! - f i) = a(Q) - a (r x - f x ) . (4.32) 

The term ct(0) corresponds to a contact term for which 
both vertices are linked to the same quark of the q-q- 
system. In intermediate steps of our calculation, a is a 
useful bookkeeping device to evaluate the diagrams in 
Fig. |l0|a. Our final result, however, will depend only on 
the dipole cross section. 

1. Inelastic one-fold differential cross section 



The diagrammatic book-keeping of the inelastic con- 
tributions to the one- fold differential cross section ( 4.3C ) 
are involved. In appendix |c[ we classify all contributing 
diagrams and we use the identities Fig. and b to show 



that many of these diagrams cancel each other. As a 
result of this analysis, the remaining non-vanishing con- 



tributions can be represented by the diagrams in Fig. 10 1 



The colour trace for the iV-th order terms of all these di- 
agrams is N c Cp\ the inelastic part of the cross section 
( 4.30[ ) is thus free of colour interference terms, it is colour 
trivial. 



(a) (l) 



(2) 



(3) 





def 

~ arbitrary number of real or contact interactions 
clef 



-e- 



+ e— e — +. 



(b) 



(c) 



FIG. 10. (a) The non-vanishing contributions to the in- 
elastic p art o f the one-fold differential photoabsorption cross 
section ( 4.3C ) as classified in appendix^, (b) Definition of 
diagrammatic shorthands used in (a), (c) Order N = 1 con- 
tribution to (al). 




We start with the N 
Fig. filial. 



1 contribution to the diagram 



L n(£) 



(I 



ri — ri — x x 



(4.33) 



This contribution is shown explicitly in Fig. |lO|c. The 
relative transverse separation ri — f\ between the quark 
in the amplitude and complex conjugated amplitude does 
not change with longitudinal position £. This is a conse- 
quence of evolving on both sides of the cut with Green's 
functions of same energy p\ . It makes the iteration of the 
N = 1 contribution to arbitrary high orders particularly 
easy. The total contribution to Fig. 0al reads 



(4.34) 



The first term is just the exponentiation of the N = 
1- contribution. The second term subtracts those TV-th 
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order contributions which involve only contact terms and 
thus do not add to the inelastic contribution. 

To evaluate the three contributions in Fig. |i"o|a2 and 
a3, we assume that the real momentum transfer occurs 
at an arbitrary but fixed longitudinal position £. The 
evolution of the free incoming F ock st ate up to £ is then 
described by the path-integral ( 4.18 ). The part of the 
expression for £ > £ is given by 



e ipM>i(l)-ri(6) e -J( d «- 

X [o- (r 2 (|) - ri(|)) + a (r 2 (£) - ri(|)) - a (0)] . (4.35) 



2. Diffractive and total photoabsorption cross section with 
one jet resolved 

The diffractive one-fold differential photoabsorption 
cross section involves by definition only contact terms. 
Due to energy momentum conservation, at least one con- 
tact term is required on the amplitude level. An arbitrary 
non-vanishing number of contact terms in the amplitude 
translates into a factor S — JC, see the discussion of Fig. |^ 
(d). Hence, the diffractive contribution to the cross sec- 
tion reads 



da 1 ^ qq n 



Here, the real interactions at £ combine to the second dadpj^ (27r) 2 



line of (4.35) and the arbitrary number of contact terms 
at £ > £ leads to the exponential of cr(0). Adding up 
all contributions of Fig. [To|a, we fin d for the differential 
photoabsorption cross section ( |4.30 ): 



dbdb<£>(h;h;a) / drdre ip 



■(r-r) 



da dp^ 



db rfb$(b;b;a) 



(2tt) 2 

*p^-(b-b) (- J L di?^l*(h-h) 



di^p- I drdf/C(b,0;r,e|A*) 



x{<r(f)+(7(r)-a(0)}/C(r,|;b,0|/i) 



x e J e 



-/^^ S (°) eI p^-(r-f) 



(4.36) 



The first term in the wide br acket s stems from the con- 
tribution Fig. [ic|al) given in ( [4.34 ), the second term de- 
notes the diagrams Fig. |Io|a2) and a3). Here, the Green's 
functions K, describe the dynamical evolution of the q-q 
dipoles up to the real interaction at £, from which point 
onwards the dynamics is given by (4.35). In analogy to 
(4.23), the lower boundary of the intergration over £ is 
determined by the last interaction point before the real 
interaction. To remove this indirectly defined variable 
from our solution, we use again the opacity expansion 
([4.20D to show that 



da dp^ 



(27T) 



db db$(b;b;a) 



dr dr 



e ip ^ r - r >jqb,0;f,£|M)/C(r,L;b,0|/i)l . (4.37) 

It is easy to check that the p^-integration of this expres- 
sion coincides with the inclusive inelastic photoabsorp- 
tion cross section ( 4.24 ). 



(5 (2) (b-r)-/C(b,0;f,L|u) 
(5 (2) (b-r) -K(r,L;b,0\n) 



(4.38) 



The one-fold differential tota l cro ss sect ion i s given by 
the sum of the contributions ( 4.37 ) and ( 4.38 ), 



da 1 ^ qq ~ 

""total 

da dpj^ 



Op 



db db$(b;b;a) 



(27T) 2 . 

ip^-(b-b) ^ e -/ £ «IC^<r(b-b) + 1 



(4.39) 



3. Non-trivial colour interference in one-fold and two-fold 
differential photoabsorption cross sections 

Colour triviality of differential cross sections cannot be 
taken for granted. It has to be established by the dia- 
grammatic techniques used in the above subsections. To 
illustrate this point, we give in the following simple ex- 
amples of untagged differential cross sections which show 
colour interference effects. We calculate the inelastic part 
of the N = 2 photoabsorption cross section (3.9) in the 
v —> oo limit in which the q-q state propagates along 
straight lines through the nuclear target. For a target of 
thickness L and homogeneous density n(£) = no, 



dn 1 ^ qq 

"°incl. 



(N = 2) 



da dp^ dp^ 

db\ db2 dbi db 2 $(bi — b 2 ; bi — S2; a) 



(27T) 



Xf 



ipj -(bi- 



b i) g _i P2 -(b2-b 2 ) 



\L 2 1 



2C 2 F 4 



[a-(bi - bi) + a(b 2 - b 2 ) 
-a(bi - b 2 ) - a(b 2 - bi)] 

Tr[aabb] B {1) + Tr[abab] B {2) 



(4.40) 
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where we have introduced the notational shorthands 

— <j(bi — bi) + a(b 2 — b 2 ) 

+a(b 1 - b 2 ) + a(b 2 - b x ) - 4 ct(0) , (4.41) 
= a(b x - b 2 ) + a(bi - b 2 ) 

-ff(bi - b 2 ) - cr(b 2 - bi) . (4.42) 



For the cross section ( 4.4C ), colour triviality is the con- 
dition that vanishes. Ho weve r, B^ ^ 0, and colour 
interference terms remain in (4.40). In contrast, the elas- 
tic part of the N = 2 photoabsorption cross section is 
colour trivial since it involves by construction only con- 
tact terms. Colour triviality (non-triviality) of the inelas- 
tic part of the photoabsorption cross section thus implies 
automatically the colour triviality (non-triviality) of the 
corresponding total cross section. 

In general, the more inclusive the cross section, the 
more likely it is colour trivial. In this sense, tagged 
cross sections are very exclusive since they require de- 
tailed knowledge of the final state of the medium. On 
the other hand, calculating more inclusive observables 



from (4.40), one recovers colour trivial examples. Espe- 



cially, the p 2 ^-integral (or pi^-integral) of (4.40) turns 



out to be colour trivial. This is ex pected, of course, since 
these int egrat ed versions of ( 4.40 ) determine the N — 2 
term of ( 4.37 ) in the v — > oo-limit. 

Based on these results, one may ask whether all one- 
fold differential photoabsorption cross sections are colour 
trivial. This is not the case as can be seen e.g. from the 
inelastic parts of the one-fold differential cross sections 



da 1 ~^ qq 

atT inol. 



dad(pf - p£) 



da- 



inel. 



dad (pi + P2) 



(4.43) 
(4.44) 



For these cross sections, the arguments of a in the ex- 
pression (4.42) for B 1 - 2 ^ are constraint by bi — bi 
■ (b 2 — b 2 ). It is easy to check that B^ ^ 0, i.e., both 
results are not colour trivial for N = 2 and thus they 
cannot be colour trivial for N = ar bitr ary. For a colour 
trivial expression, it is crucial that (3J3) is not integrated 
over some average or relative momentum but over the 
transverse momentum of a final state particle. This is 
necessary for using the identities of Fig. I (a) and (b). 



F. Mueller's dipole formulas as eikonal 
approximation 

Explicit calculations of the pj_-integrated and one-fold 
differ entia l photo absorption cross sections ( 1.24 )-( [!. 26| ) 
and (4.37)- (4.40) respectively require an explicit repre- 
sentation of the path integral K,. This necessarily involves 
an approximation. In this section, we study the eikonal 



approximation which substitutes the path integrals K, by 
their ultrarelativistic limit 



lim K(L, r; 0, b\(i) = <S (2) (b - r) 5(b) . 



(4.45) 



We consider a homogeneous density distribution n(£) = 
no for a target of longitudinal extension L and we write 
for the Glauber-type suppression factor 



5(b) 



:L<r(h) 



(4.46) 



Also, we neglect the phase factors exp{iq (£ — £')} in this 
v — > oo-limit, see our discussion of ( gjf ), flOg ) above. 
In the ultrarelativistic limit, all photoabsorption cross 
sections reduce to standard Glauber- type expressions: 



da J db $ (b; b; a) 



[1 - 5 2 (b)] 



(4.47) 



J da J db$(b;b;a) 



[1-S(b)f 



(4.48) 



lim a 1 ~* qq 



2N c a cm J da J db$(b;b;a) 



[l-5(b)] 



(4.49) 



This total photoabsorption cross section is Mueller's 
dipole formula first derived in Ref. [[5). Also, the elastic 
cross section is known to describe the diffractive contri- 
bution to the deep inelastic scattering structure function 
F2 in the ultrarelativistic limit [ [| . 

The one-fold differential cross sections read in the 
eikonal limit 

lim da ^ q " = f dbdb<f>(b;b;ayrt< h -V 

f-><» da dpf (2tt) J 

x [S(b - b) - 5(b) 5(b)] , 

da2;»T qq a em N, 



lim ~ diff - 1 - 

f->oo da dpf (2tt) 



(4.50) 

db db$(b;b;a)e 4p ^ (b - 5) 



[1 - 5(b)] [1 - 5(b)] 



(4.51) 



dal^^a^N, 1 dbrf b$(b;b;a)e^ b bi 
y^ca da dpf (27r) 



lim 



x 



[l + 5(b b) 5(b) - 5(b)] . (4.52) 



The elastic contribution ( 4.51 ) has been obtained pre- 
viously in calculations of the diffractive component of 
DIS electron nucleon scattering with one resolved jet in 
the final state [Hi. Furthermore, the total cross section 



( 1.52 ) was previously obtained by Mueller (see eq. (27) 
of Ref. [ |8|) from a one-loop calculation arguing then 
for its general validity. In passing, we r ecall Mueller's 
interpretation of the four terms in ( 4.52 ): for the case 
Pi *C Q 2 , the first two terms turn out to be of equal 
size while the last two terms are negligible [ pf . Looking 
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at the diagrammatic contributions one concludes that in 
this limit the term proportional S(h — b) gives the prob- 
ability that a quark which gets many random kicks is 
found with relatively small transverse momentum. The 
term proportional to 1, on the other hand, corresponds 
to the case of no scattering and can be viewed as the 
quantum mechanical shadow of the term proportional to 
5(b-b)[|. 



G. Corrections to the eikonal approximation 

To go beyond the eikonal limit, we discuss now the 
Gaussian dipole approximation for the path-integrals /C. 
It is based on the observation that the main support in 



the path integral (4.19) comes from small transverse dis- 
tances r = |r|, where the cross section <r(r) in (4.8) has 
a leading quadratic dependence: 



o-(r) w C(r)r 2 



(4.53) 



In explicit model calculations, one finds [ [H| that the 
r-dependence of C(r) is a slow logarithmic one, and can 
be neglected. For the dipole cross section in (4.8), an 
expansion to order r 2 confirms this feature. One finds 



|a (/«±)| 2 , 



(4.54) 



where the Kj_-integral depends logarithmically on the ul- 
traviolet cut-off k c which one has to introduce in this 
approximation. We note that C provides a measure 
of the average transverse momentum transfer (Kj_) in 
a single scattering. For sufficiently small r, when the 
r-dependence of C can be neglected, the path integral 
(4.19) reduces to that of a harmonic oscillator [ |i~5| ] 



^4 

/C osz (r 2 ,£;ri,0|/^) = — exp {iAB{v\ + r 

-2zAri • r 2 } , 



2 sin(fi L) ' 
B = cos(^L) , 



with the oscillator frequency 



n = 



1 — i hqC 



^/2 



(4.55) 
(4.56) 
(4.57) 

(4.58) 



For what follows, it is important that the (i — > oo-limit 
of /C osz coincides with the eikonal limit ( 4.45[) of the un- 
approximated path-integral. To see this, we expand the 
norm of JC osz to leading order in /i and the phase to next 
to leading order, 



/C osz (r 2 ,£;ri,0|/z) 



m + o(m°) 

2iriL 



LtIqC r o n 

— rf + r z 



+ r e -r] +0( M " 1 ) 



6& (r e -r) e-^» iC - 2 , 



(4.59) 



which coincides with (4.45) for the quadratic ansatz 
a (r) = Cr 2 . The harmonic oscillator approximation 
thus provides an explicit representation for the path- 
integral JC which preserves the correct high energy limit. 
At least numerically, this allows for an explicit study of 
the 1 / /x-corrections to the eikonal limit of the cross sec- 
tions ( |4.24| )-(p~2e|) and fl4.37| )- fl4.39| ). To explore some 
qualitative features of these l//Lt-corrections analytically, 
we take here recourse to a Gaussian approximation $g 
of the incoming Born wavefunction: 



$ G (r,r;a) = ip(r)i/j*(r) , 



V.(r) 
1 



ttK- 



■ exp (—r 2 /R 2 



— = (1 - a )aQ 2 



(4.60) 
(4.61) 

(4.62) 



where the radius R is chosen to reproduce the character- 
istic widt h of the f all- of f of the exact Born probability $ 
given in (B8) and (B14). 

The time-evolved final state wavefunction reads 



*/(r) = J dr/C osz (r,L;r,0|^)V(r) 



ttR 2 



■ exp 



r 

R? 

r.2 



1 



n LCf 2 



i—c 1 -—c 2 + 0{\i 6 ) 



/./, 



iqL 
(4.63) 



Here, we have used the phase convention ( 4.27 ) to obtain 
the term iqL in the exponent. The shorthands c\ and 
C2 denote the leading real and imaginary l//Lt-corrections 
to the final state wavefunction: 



Cl= 6 



C2= 15 



n 2 L 2 C' z + 6 



n LC 



12 

R* 



n 3 L 3 C 3 



10- 



]L 2 C 2 



40 



n LC 



60 
~R? 



(4.64) 
(4.65) 



Deviations from the eikonal limit can be calculated with 
the help of ( |4.63 ). For example, we find for the total 
photoabsorption cross section (4.26): 



1 ^99 _ l: 7 "^99 

"total J^ 1 ^ "total 



J da J dr$ G (r, r; a) S(r) 



x2 



1 



e" r C2/ ^ cos 



-Ci + qL 



(4.66) 



In the v — ► cxi-limit, this difference vanishes by construc- 
tion. To understand which scales determine the devia- 
tions from the eikonal limit, we recall that l/(/ii? 2 ) = q 
has an interpretation as inverse coherence length 1/lf 
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and uqLC is the total transverse momentum squared ac- 
cumulated during the rescattering over a distance L. The 
eikonal approximation is thus justified if the phases de- 
termined by the coherence length and the total transverse 
energy E^ 1 are negligible: 



qL < 1, 

E tot L= ^LC L ^ 1 
2/i 



(4.67) 
(4.68) 



We conclude this section with two technical remarks: 
1. The explicit form of the incoming q-q wavefunctions 



(B6), (B13) is given in terms of the Bessel function Kq 
which has an integral representation in terms of a Gaus- 
sian in r, 



Ko(e\r\) = 



1 



(4.69) 



The evolution of the Gaussian wavefunction ip(r) in 
(4.63) can thus be used for the exact calculation if l/R 2 



is taken to be the width of the integrand of ( 4. 69] ) 
and the ul -integration is done afterwards. Howev er, to 
order terms proportional to l/R 4 appear in ( 4.63j ) 
and thus the tjl -integration cannot be done analytically 
if l//i-corrections are included. Only the /i — > oo-limit 
can be accessed analytically in this way. 

2. A surprising difficulty occurs in the calculation of 
the photoabsorption cross secti on (4.24 ) from the har- 
monic oscillator approximation (4.55) if one tries to do 
the r e -integration first: 

F osz (r, r) = J dr e /C OS z(r, 0;r e ,L|^) /C osz (r e , L; r, 0|/i) . 

(4.70) 

After a lengthy but straightforward calculation, we find 
for the fj, — > oo-limit 

lim F osz {f,r) = 8 {2) (f-r) e-i noLCr2 , (4.71) 

fl — >oo 



whose exponent differs from that of the eikonal dipole 
formula (4.45) by a factor i. This mismatch is an arte- 
fact stemming from a calculation which performs sim- 
plifications on the cross section level before completing 
the dynamical evolution on the amplitude level. Start- 
ing from ( 4.63 ) for \&/(r e ) an d calc ulating \P/(r e ) ^fj(r e ), 
one obviously reproduces for (4.24) the eikonal expression 
(4.47) in the fx — > oo-limit. 



V. PHOTODISSOCIATION VIA NON-ABELIAN 
STOKES'S THEOREM 

In this section, we argue that the dipole cross section 
(4.8) parametrizes the transverse components of the chro- 
moelectric field strength correlations (F F) in the nuclear 



medium. This suggests an interpretation of c(r) which 
holds model-independent, i.e., irrespective of whether we 
have build up ( [Q| ) from the model ansatz ( 2.10Q , or from 
some other model parametrization of the colour target 
field. 

We first recall the non-abelian Stokes's theorem. This 
relates the integral of the non-abelian vector potential 
along a closed contour C to the area integral of the 



field strength tensor [ 
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TtV exp \ J> dz^A^iz) 

= TiV alea exp ( f da^iy) U Ay F^(y) U yA ) . (5.1) 



Here, the Ua v are parallel transporters. They ensure 
that the r.h.s. is gauge-invariant. Equation (5.1) does 
not depend on the choice of the position A, since cyclic- 
ity of the trace allows to change the reference point by 
substituting U A y — * Uba UAy The area ordering Parea 
is defined by disecting the area enclosed by C in many 
(ultimately: infinitcssimally small) areas in such an order 
that the parallel transporters C% along the infinitessimal 
areas combine to C. As a consequence of the non-Abelian 
Gauss theorem, (5.1) turns out to be independent of the 
orientation of the surface da IJtv {y) [ [361. 



The total photoabsorption cross section (4.1) is closely 
related to the Stokes's theorem, as can be seen from the 
representation 



7 lt7i qq = "cm J da J dbi db 2 dbi dh 2 $(Ab; Ab; a) 



TrV exp 



' C(r 1 ,r 2 ,ri,r 2 ) 



dzpA»(z) ) (5.2) 



Here, C is defined by the transverse paths i"i(£), Fi(£), 
r 2(£)> ^a(0 which denote the positions of the rescattering 
quark and anti-quark in the path integral representation 
(|2.17| ) of the corresponding Green's functions. We can 
consider C as a closed path since we work for a static 
potential A^ — 5^ A for which the Wilson line along 
purely transverse directions dxj_ A±(x) vanishes. The 
photoabsorption cross section (^) is thus given by a 
closed Wilson loop which sums up the gluon field strength 
of the medium in the area determined by the paths of the 
rescattering q- and q- quarks. 

The connection between total hadronic cross sections 
as fl5.2| ) and closed Wilson loops has been studied exten- 
sively, see e.g. [ 26, Phenomenological applications 



of the non-abelian Stokes's theorem (5.1) were pioneered 
by Dosch and coauthors [ B% |38] . In Dosch's approach, 
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one takes recourse to an axial gauge in which the par- 
allel transporters on the r.h.s. of fl5.1| ) reduce to unit 
operators. The remaining express ion is then expanded in 
powers of F. Our discussion of (|5.2|) borrows from this 
strategy. Clearly, we cannot choose an axial gauge to re- 
move the parallel transporters in (5.1), since our gauge 
freedom is already exhausted by the definition ( 2.10 ). 
However, we can invoke the proof of colour triviality of 
( |5.2| ) to discuss the equivalent abelian problem with a 
rescaled vector potential A, t — > \JCf A^ . We find 



ds 



exp I— i g \/Cf J> dsA^(uj s 

I ( ■ fn~ f e. r \ duJ s duJ s , . 
= ^exp\ig^C F J F tlu {Lu s )-g--Q—dsdx a 

r z 2 /t 2 (S) 

-i g y/C^ / d£ / dx\ El (x ± , £) 
Jzi J ri (0 



= 1 - 



ri(0 



x J B^(x x ,0^ ± (x'x,e + 0}+O(n 2 (0) • (5.3) 

Here, ui£ denotes the path around the contour C, and 
the notation is specified in Fig. [H]. 




r 2 (z) 



FIG. 11. The closed Wilson loo p b etween longitudinal 
positions z\ and Z2, as considered in (5.3). 



We have used the particular form of a static colour po- 
tential A^ = 8^0 Aq to rewrite the r.h.s. of (5.1) in terms 



of transverse electric fields. The detailed arguments for 
the shift of integration variables needed to arrive at the 
last equation of ( |5.3D are given in appendix B of Ref. [ 
E2| . The same appendix contains the derivation of the 



dipole cross section (4.8) from a closed abelian Wilson 
loop: 

cxp f-i g VCV j> ds i^Us) ^ i 
cxp \-igy/C^£ <% [j4o(ri(O,0-4>(ra(O,O] 



exp 



(5.4) 



section to the two-point correlation function of the trans- 
verse electric field strength, 

<7(ri(0-r 2 (0) = ( d£ dx\dx'l 

x^( X± ,e)^(x' ± ,e + o). (5.5) 



As long as we parametrize this dipole cross section by the 
Gaussian approximation <r(r) = Cr 2 , only the average of 
the colour field strength enters our final results in form 
of one single parameter C. The final result is model- 
independent in the sense that it does not depend on the 
model-specific way in which the colour field strength giv- 
ing rise to C was modelled. 



VI. CONCLUSION 

Colour triviality provides a crucial simplification of an 
otherwise untractable problem. In general, the iV-fold 
rescattering of a parton in a colour target field leads in 
the cross section to colour traces over 2N generators. For 
sufficiently exclusive processes, the huge number of colour 
interference terms thus limits explicit calculations of the 
soft nuclear dependence of hard partonic processes to the 
case of very few rescatterings [ M. Colour triviality is 
the consequence of a complete diagrammatic cancellation 
between these different colour interference terms. For 
colour trivial observables, all non-abelian complications 
reduce to the rescaling of coupling constants by appropri- 
ate powers of Casimirs. The multiple scattering problem 
thus becomes an abelian one. The soft nuclear depen- 
dence of hard colour-trivial observables can be described 
in terms of a QCD dipole cross section which absorbs the 
leading medium-dependence in a one-parameter estimate 
C of the average target colour field strength. This makes 
the identification of colour-trivial observables of particu- 
lar interest for relativistic heavy ion collisions at RHIC 
and LHC where very little is known about the medium 
a priori but many "hard probes" are expected to receive 
sizeable nuclear modifications due to rescattering effects. 

In the present work, we have derived an explicit general 
expression for the rescattering effects on a hard coloured 
parton inside a spatially extended nuclear medium. We 
have then studied in detail for the simple example of vir- 
tual photodissociation under which conditions this gen- 
eral rescattering formula results in colour-trivial observ- 
ables which can be parametrized by a QCD dipole cross 
section. In detail: 



Comparing (5.4) to (|5.3|), we can relate the dipole cross 



The non-abelian Furry wavefunction ( 2. IS ) derived in 
section || is a high-energy approximation to the solution 
of the non-abelian Dirac equation in a spatially extended 
colour field. In contrast to other approximate solutions [ 
it is accurate up to order 0(l/E 2 ) in the phase. This 
is known to be indespensable for calculating the nuclear 
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dependence of observables which are determined by the 
destructive interference between different production am- 
plitudes as e.g. the LPM-effect or the nuclear d epen dence 
of Drell-Yan yields. The Furry wavefunction ( 2.1 2| ) thus 
provides a unified starting point for the description of the 
nuclear dependence of a large class of observables. 

Based on the Green's function G which describes the 
dynamical evolution of the non-abelian Furry wavefunc- 
tion, we have derived in section IV a set of diagrammatic 
identities which play the key role in proofs of colour triv- 
iality. These identities exploit the integration over the 
transverse momentum of a single final state particle only. 
They are thus much stronger than diagrammatic cancel- 
lations implied by the optical theorem. As a consequence, 
they ensure for the example of the 7* — ► qq process a 
colour trivial result not only for the total, elastic and 
inelastic inclusive photoabsorption cross section but also 
for cases in which one jet is resolved in the final state. All 
these observables can be described in terms of the same 
one-parameter QCD dipole cross section. 

Earlier applications of these diagrammatic techniques 
as well as statements about the colour triviality of the 



photoabsorption cross section exist 

H |J. Here, we have given relatively short and complete 
proofs of these statements by exploiting the advantages of 
a new and compact configuration-space notation implied 
by the Furry wavefunction. In contrast to previous dis- 
cussions, our formulation includes the transverse dynam- 
ical evolution of the partons. For the example of the pho- 
toabsorption cross section studied here, this allowed us 
to quantify the leading deviations from the eikonal limit 
which grow at fixed energy proportional to L 2 . More- 
over, the compactness of our notation made it possible 
to discuss in detail technically rather involved processes 
as e.g. the inelastic and total photoabsorption cross sec- 
tion with one jet resolved in the final state. This points 
to the strength of the present approach which we believe 
to be suited for the explicit discussion of more compli- 
cated processes as, e.g., further studies of the transverse 
momentum dependence of the non-abelian LPM-effect [ 
|l8| , |22| ] or photodissociation including initial state gluon 
radiation [ p3| . 
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APPENDIX A: NON-ABELIAN FURRY 
APPROXIMATION 

In this app endix , we derive the non-abelian Furry 
wavefunction ( 2.12 ) from the set of Af-scattering Feyn- 
man diagrams ( 2.11 ). For an Y-fold rescattering dia- 
gram, we use the notation Piv+i = P for th e final state 
momentum. To simplify I^ N \y) in ( 2.11 ), we do the 
longitudinal momentum integrals by contour integration 



(27r) p 2 — m 2 + ie 
_ (fa + to) 70 



2 p\ 



-1) 



(Al) 



We note that, strictly speaking, the position variables Xi 
in this expression are integration variables and do not 
coincide with the center of the i-th scattering potential. 
A more detailed analysis of ( 2.11 ) involves for a particu- 
lar model also the study of the pf -poles of the Fourier- 
transformed single scattering potentials. For a Yukawa- 
type potential l/[(pi — Pi-1) 2 + M 2 ], e.g., these poles 
give additional contributions to the -integrals, which 
are however exponentially suppressed due to the Debye- 
screening mass M. These details are discussed explicitly 
in Refs. [ |l3|, [l6| and for contact terms in Ref. [ ^] . In the 
end, the only 0(l/i?)-contribution to the -integration 



turns out to be given by (Al). On the r.h.s. of this equa- 



tion, pf is determined by the pole value to order 0(1/E), 
p\ = p — t| — where p 2 = E 2 — to 2 . Using xq = y, we 



can rewrite (2.11) as 



/ n r ,2 ± 
- (A + m)70 H^(x,:)]e 4 ^-^- x -) 



2 pf 



« w (p). (A2) 



We now consider the spinor structure of this expression. 
To order 0(1/ E 2 ), each quark propagator introduces a 
numerator 



(J) t . + to)7q ps £(70 - 73)70 - 7^ • pf 70 + 73 7o 



2p 



(A3) 



The normalization of I^ N \y) is needed to leading order 
in energy only. This allows us to neglect the mass term in 
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( |A3| ) . For the same reason, it is sufficient to keep only the 
leading term £"(70 — 73)70 for all numerators (A3) with 
i > 2. Further simplification of (A2) is then possible, 
using 

(£(7o - 73)70)" = 2"- 1 E n ( 7o - 73)70 . (A4) 

The case i = 1 is different: depending on the explicit 
form of the production vertex, the leading order contri- 
bution of this numerator can cancel. We have to keep the 
numerator ( /p \ + m)7o to order 0(1/ E). To this aim, we 
substitute in (|A2| ) 

d 



(j6i+m)7o — ► £(70 - 73)70 - h ■ g-Jo , (A5) 
where the differential operator acts on the big bracket 



in (A2). Equation (A5) is equivalent to (A3) with the 
transverse components written in configuration space. In 
a coordinate system with p || n, the above equation takes 
the form 



hlo = p(lo ~ 73)70 + a 1 " • (pi 



-a 



2p 2p 



— > P (7o - 73)70 +icn ■ a-(p-pn). (A6) 

oy 

In the same coordinate system, 0.5 (1— 7370) v(p) = v(p). 
Acting with (A6) on v(p), we thus find to leading order 
in E an expression proportional to Dv(p). With these 
steps, the amplitude (A2) takes the form 



I {N) (y) = e-' ipyL DV 



'i-lJ 



Y[[-iA (xi)]Go(xi-i,Xi\p)\ J dx A 



x Go (xjy , xlpjfoo (xi , x L , p) (p) . ( A7) 

Here, we have introduced two elements of the abelian 
Furry approximation: the outgoing transverse plane wave 
Foe and the free Green's function Go of ( |2.15 ) which we 
have identified here with the Gaussian p ^--integrals, 

G (x^ 1 ,a;f_ 1 ;x^,a;f \p) 

(2tt) 2 ' [ ' 



For the path-ordered product in (A7) which involves 
these Green's functions, we introduce the shorthand 



G (JV) (y±,yL;xj_,x L |p) 



x i-i) \ Go(y±,yL;xi,x±\p) 



V I f[ Jd 3 Xl Q(x 

\ i=l / 

f IlH M*i)] Go(xl,xf; X j- +1 ,xf +1 \p) ^ 



(A9) 



where xjy+i = x^. Allowing for arbitrary many gluon 
exchanges, we have to sum over N. 



G(y±,y L ;^±,x L \p) = ^2 G {N) (y ± ,y L ;x ± ,x L \p) . 

(A10) 



N=0 



From this one sees easily that G^ N ' corresponds exactly 
to the N-th order O(Ag) term in ( |2.14| ). The sum ( ]2.12| ) 
over iV-fold scattering diagrams takes the form of a non- 
abelian extension of the Furry approximation. 

We emphasize that the approximations used in our 
derivation and especially the different treatment of the 
quark propagators at the production vertex and in the 
final state rescattering part of the amplitude 1^ leads 
to a consistent high energy expansion with a norm accu- 
rate to leading order in 0(1/ E) and a phase factor ac- 
curate up to order 0(1/E 2 ). The same approximations 
were employed in recent calculations of rescattering am- 
plitudes for which the production vertex P is a photon 
or gluon emission vertex [ Ol p7| . 



APPENDIX B: QQ FOCK STATES 

In this appendix, we give details of the derivation of 
transverse and longitudinal components of the squared 
incoming Fock state 

$(Ab; Ab; a) = $±(Ab; Ab; a) + $ L (Ab; Ab; a) . 

(Bl) 

We start with the longitudinal polarization , satisfying 
e L ■ k = and e l2 = 1, 



1 



(B2) 



The corresponding vertex function reads to leading order 
in energy 



r = e 



' y/(l — a)ar 8 r y 



(B3) 



where r, r' are the helicities of the quark and antiquark. 
The square of this longitudinal vertex function, summed 
over the spin of the final state particles, reads 

^r I r i * = 2Q 2 (l-a)a. (B4) 

r.r' 

This is a kinematical prefactor which factorizes in the in- 



tegrand of (3.7). For the firs t interaction in the photodis- 



sociation amplitude in (3.7), we t ake some longitudinal 
position z a . The y-integration in (3/7) can then be done 
analytically. We find 
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I L (b 2 - bi\y L ,z a ) 

d 2 yj_ G (b 2 ,z a ; y\p 2 ) G (y; bi, 2 a |Pi) 



/' 



2iri (z a - y L ) 



exp • 



iH (b 2 - bi) 

• 2{z a -y L ) 



■}, (B5) 



where jtx = (1 — a)a The virtual photon has to dissoci- 
ate at longitudinal position < z a in order to make an 
interaction at z a possible. This limits the range of the 
t/i-integral: 

Z a 

I L (Ab\z a )= J dy L I L (b 2 - bi|y L ,z a ) 

— oo 

(1 — a) av 



2-ki 



2K (e|Ab| 



e= V(l-a)«<3 2 



(B6) 
(B7) 



If the first interaction occurs at z a in both the amplitude 
and complex conjugate amplitude, then the square of the 
incoming longitudinal qq Fock wavefunction is defined as 



the combination of the squared emission vertex (B4) and 
the free time evolution (B5), 



$ L (Ab; Ab;a) 



4j/ 2 a(l 
Q 2 a 2 (l-a) 2 



1 ' ' [ ^ 7 L (Ab|z Q )/£(Ab|z a ) 



2 (2tt) 2 



4X (e|Ab|)X (e|Ab|). (B8) 



This expression appears for longitudinal polarization in 
the photodissociation cross section 



For our calculations in section IV, we have to consider 



the more general case that the first interaction vertex 
occurs at z a in the complex conjugated amplitude Mg, 
but at z b > z a in the amplitude Mg. This implies the 
further free evolution of (JBq) from z a to z b , 



I L (b 2 - bi\y L , z b ) 

= J d 2 b\ d 2 b' 2 Go{b 2l z b -b' 2l z a \p2) 

X I L (b' 2 - b[ | y L , z a ) G (b[ , z a ; b x , Zb\pi ) 
= J dr b (z a ) K (r b (z b ) 7 z b ; r b (z a ), z a \va{\ - a)) 

xI L (r b {z a )\y L ,z a ) . (B9) 

Here, we have used the notation introduced in equations 
( [hl6D - ( [4.19D with r b (z a ) = b' 2 - b[ and r b (z b ) = b 2 - 
b\. Equation ( p39| ) allows us to circumvent a notational 
problem which stems from the identity: 



dy L I(Ab\y L ,z b ) 



a iq(z b -z a ) 



dy L I{Ab\y L , 



(BIO) 



Equation (BIO) shows that if the points of first in- 
teraction differ in AfJ and Mg, then phase factors 
exp (i q (z b — z a )) arise in the squared free incoming wave- 
function ( p38| ). Hence, strictly speaking, we cannot do 
the j/L-integral before specifying the first points of in- 
teraction. Previous discussions of the q-q-dipole do not 
have this difficulty, since they neglect these phase fac- 
tors - a praxis which is justified for nuclear targets 
of size L <C 1/q. If we want to keep the phases 
exp (i q ( z b — z a )) and yet write the photoabsorption cross 
section ( |3.9|) in terms of we can use ( |B9| ) for a 
simple convention: whenever the free Green's function 
Kq (r b (z b ), z b ;r b (z a ), z a \va(l — a)) acts on the first ar- 
gument Ab of $i(Ab;Ab;a), the result is a phase 
exp{iq(z b — z a )), whenever it acts on the second argu- 
ment, the result is the complex conjugated phase. This 
allows us to expand the cross section ( 4.23 ) in powers of 
the opacity without neglecting phase factors and without 
giving up the simple representation of cr 7 in terms 
of 

We now turn to the transverse polarizations 



c£(A) = — (0,MA,0) . 



1 

72 

rI(y) = e^(A)r" 

-1 



2va 



(Bll) 

[5 r yiV y ■ e ± (X) (r'(l - 2a) + A) 



+—5_ rtrl (l + \r')} 



(B12) 



In contrast to the longitudinal case, this vertex does 
not factorize in the integrand of (3.7). However, the y- 
integration can be done in analogy to (B6), 



^(Ab|z a ) 



6 a 

= J dy L |rf 2 yie i?!,l " £|yi| Go(b2;y|p2) 



xrI , (y)G (y;b 1 | Pl ) 
(1 - a)av ^ lQZa ^ rT 



2-ki 



2ri(Ab)^ (e|Ab|). (B13) 



The squared transverse wavefunction reads then 

2 I/A,r,r' 



$ ± (Ab; Ab;a) = 



Av 2 a(l — a) 



l£(Ab\z a )l£ (Ab\z a ) 



Ab - Ab 



2(2tt) 2 |Ab| |Ab| 



KxieAb) ifi(eAb) (a 2 + (1 - a) 2 ) 



2 (2tt) ; 



■ Ko(lAb) K {eAb) 



(B14) 



For the case that the points of first interaction z a and 
z b are different in the amplitude and complex conjugated 
amplitude, the convention explained below (B1C) applies. 

The arguments of th e modified Bessel functions K$ and 
K\ in (B8) and (B14) specify the transverse separation 
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of the q-q dipole pair. This separation is given approxi- 
mately by 1/e = 1/ \Z(T — a) a Q 2 + m 2 . The higher the 
virtuality of the photon, the smaller the transverse size of 
this dipole. In general, the approximation of the virtual 
photon by the q-q Fock state is reasonable at sufficiently 
high virtuality, Q 2 3> 1 GeV 2 , above the vector meson 
dominance region [[7j. There, the q-q separation is small 
enough to render the effects from string tension negligible 
and the quarks undergo independent scatterings. For a 
recent attempt to model interactions between the quarks, 
see Ref. \ M. 



APPENDIX C: CLASSIFICATION OF 



DIAGRAMS FOR (4.30) 



In this appendix, we use the identities Fig. and b 
to simplify the diagrammatic contributions to the inelas- 
tic part of the differential photoabsorption cross section 
(4.30). For each iV-th order term in the opacity expan- 



sion of (4.30), we use the notational shorthand 



N 



IK 

i=l 



(CI) 



Here, the index i labels the scattering potentials linked 
to the N-th order term, i = 1 being the first scattering 
potential. The superscripts m, mi, o denote how the i-th 
scattering center is connected to the g-g-system. Each 
scattering center is linked twice to the g-g-system, once 
a momentum k±_ flows into the system, once it flows 
out. rii € [1,2] specifies whether the inflowing momen- 
tum is transfered to the upper quark line of momentum 
Pi {ni — 1) or the lower quark line of momentum p2 
(n-i = 2). Analogously, m t G [1,2] specifies whether the 
outflowing momentum comes from the upper (to; = 1) 
or lower (m^ = 2) quark line. The remaining superscript 
o G [r, v, w] specifies the position of the two vertices from 
the i-th interaction w.r.t. the cut: either both vertices 
stand to the right of the cut (o = v) or to the left of the 
cut (o = w), or they are a real contribution (o = r) where 
one vertex stands to the right and the other to the left 
of the cut. We illustrate this notation with the examples 
in Fig. HI 



A',"A7 = 



. ll.w . 12,r 

A, A 2 = 




FIG. 12. Example of our notation (Gl) for diagrammatic 
contributions to the photoabsorption cross section. For more 
details, see the text following (|Cl|), 



All iV-th order contributions to the differential pho- 



toabsorption cross section (4.30) fall into exactly one of 
the following three classes: 

1. There is at least one factor A- ' T in the diagram 



(El) 



2. There is at least one factor A l2 ' r or A- ,r in the 
diagram (CI), but no factor A.' r (arbritrary k) 



and no factor A 22,1 ' with k < j. 
3. (a) There is no factor A^ 1,r (arbritrary k) but at 



least one factor A 12,r or A 21,r in the diagram (CI), 



and there is at least one factor A k ' r with k < j. 



(b) There are no factors A 
diagram (CI). 



11, r 



A 12,r A 21,r ■ ,i 

A ■ or A- m the 



We now use the identities of Fig. |7]a,b to show that many 
of the diagrams in these three classes cancel each other. 
In this way we determine the only remaining contribution 
for each of the three classes: 

Class 1: The only non-vanishing contributions are 



those which contain no terms A 



A, or A l with 



o G [r, v,w\. They are shown in Fig. |lfj((al). 
Argument: Consider the subclass of diagrams containing 



factors A 



22, o 



,.12, o A 2l, o 

A or A i ■ 



In case that there is more 



than one such factor in the diagram, consider the term 
with the largest index I, I = k say. Leaving all terms 
j =/= k unchanged, we find in this class of diagrams ex- 
actly one diagram with A 22 ^ , A 22 ^ and A 22 '™ on the A:-th 
position. These three diagrams cancel each other due to 
the identity Fig. ^a. Also, for all terms j ^ k unchanged, 
we find exactly one diagram with A 2 , 1,1 ' and A k ,v , which 
cancel due to the identity Fig. [?]b. For the same reason, 
the two diagrams with A 12 ' r and A 12,w cancel each other 



|12,M) 



(note that A k 
in the same diagram) 
tion which contains terms A 22, ° 



= A k ' specifies the same k-th term 



21,to 

As a consequence, no contribu 
\ 2 - or Af 



a12,o a2\,o ■ 

A, ' or A, ' gives a 



non-vanishing contribution. 

Class 2: The only non-vanishing contributions contain 
exactly one real term A 12 ' 1 or A. ,r with no real term 

A]} ,r or A 22,r , k arbitrary, and with no contact terms 



k k 
J^2,v j^2 



A 2 , 2,v , or Alf' w for k > j. They are shown 



,22,ui 



I ^(a2) 



in Fig 

Argument: Choose in each diagram the term A X2 ' X or 
A 21 ' x with the lowest index j. Consider diagrams which 



contain terms A 



22. o 



Af' or A k 



21, o 



with k > j and o 



r,v,w. Taking k maximal and leaving all other terms 
unchanged, these diagrams cancel due to the identities 
Fig. 0a and b. 

Class 3: The only non-vanishing contributions contain 
exactly one real term A 22,r but no real terms Aj ,r , A 12,r , 



or A 21 ' r and no contact terms A 22 '' 



Al 



with k > j. They are shown in Fig. |10|(a3) 



A 12 ' v orA 12 > w 



22 



Argument: Consider in each diagram the term A 22,r with 
lowest index j. Look then at the largest index k > j 
linking to the j>2-hne. Leaving all other terms unchanged, 
the sum of the contributions with different configuration 
at position k ensures cancellation: A 2 k 
cancel each other due to identity Fig. f?]a. A 
A 12 ^ , A 21 ' w cancel each other due to identity Fig. f?]b. 
Thus only contributions with exactly one real term A - ' r 
survive cancellation. 



Af> v andAf'™ 

21, r ,21.v , 
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